Real and complex random neutrino mass matrices and 613 
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Recently it has been shown that one of the basic parameters of the neutrino sector, so called 
613 angle is very small, but quite probably non-zero. We argue that the small value of O13 can 
still be reproduced easily by a wide spectrum of randomly generated models of neutrino masses. 
For that we consider real and complex neutrino mass matrices, also including sterile neutrinos. A 
qualitative difference between results for real and complex mass matrices in the region of small ^13 
values is observed. We show that statistically the present experimental data prefers random models 
of neutrino masses with sterile neutrinos. 
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Since their first appearance in modern science neutri- 
nos have played an important part in our understanding 
of the laws of particle physics. To see how important the 
weak interaction of neutrinos is in Nature lets just men- 
tion the mechanism in which the Sun is shining [T]. It 
is also very well known that they influence evolution of 
the whole Universe as their tiny masses have an impact 
on the dynamics of the expansion of the Universe [2J [3] . 
Without any doubt the investigation of the properties of 
these particles can reveal many interesting, hidden until 
now physical phenomena or explain many hypothetical 
ideas. 

The theory of neutrino oscillations requires nonzero 
neutrino masses as well as nonzero neutrino mixing an- 
gles. 

The T2K Collaboration recently announced that a new 
measurement [4] has yielded a nonzero ^13 for CP con- 
serving case with Sep = (results at 90% confidence 
level), i.e. 

0.03 < sin^ 2(?i3 < 0.28, normal mass hierarchy, (1) 
0.04 < sin^ 2013 < 0.34, inverted mass hierarchy. (2) 

The Minos collaboration published similar results few 
days later [5]. These results with nonzero ^13 are of pri- 
mary importance for the future of particle physics [BJ [7] . 

Present global fits of the data give [S] (at 1 cr C.L., 

m) 

7.05 • 10"^ eV^ < Am^i < 8.34 • lO"'"^ eV^ (3) 

2.07 • 10-3 eV^ < Aml^ < 2.75 • 10"^ eV^ (4) 

0.39 < sin^ 023 < 0.5, (5) 

0.291 < sin^ 012 < 0.324, (6) 

0.008 < sin^ 013 < 0.036. (7) 

There are many theoretical models which tune the 
space of possible neutrino parameters to get agreement 
with these experimental data. The new data from the 
T2K collaboration generated a new challenge for that. 
Frequently the pattern of a flavour symmetry is invoked, 
for recent theoretical adjustments of the nonzero ^13 mix- 
ing angle, see e.g. [10]. Here we face the problem differ- 
ently asking how natural are small, nonzero values of ^13. 



Probability 




10-5 10-"^ 10-^^ 10-2 10-' 1 



FIG. 1: Frequency probability for R parameter obtained 
from 10^ randomly generated three dimensional different 
types of matrices. No constraints given by Eqs|3][7] are ap- 
plied. 

Naturalness means that the neutrino mass matrices are 
generated randomly. Apart from neutrino physics, see 
e.g. [Ill I17j . random matrices came to be an eminent 
tool used in many different fields of science, including 
day-to-day life problems [TH]. 

Let us first consider a symmetric, three dimensional 
matrix with random elements in the range [—1, 1]. Phys- 
ically this case realizes Majorana type of neutrinos. Simi- 
larly, we consider also Dirac neutrinos (generated by gen- 
eral real random matrices) and so-called see-saw neutri- 
nos with masses rui, defined through relation rnj^M^^mD 
(here again matrices m d and rriR are generated randomly 
in the range [—1, 1]). Next, we apply the singular value 
decomposition theorem |19| to calculate the mixing ma- 
trix U which diagonalizes the random matrices we started 
from. Using the standard parametrisation, we are able to 
link experimental mixing angles to the elements of matrix 
U, e.g. sin2 0i3 = \Ui3?. 

In Figjljwe recover one of the results discussed already 
in [TT], see Fig.l. Here R = Ato|]^/Ato§2 is the ratio of 
the smallest to the next smallest values of differences of 
squared neutrino masses, mi < TO2 < 7713. Both in Figjl] 
and next plots results are generated based on 10^ random 
mass matrices. Hall et al. approach has an advantage of 
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FIG. 2: Frequency probability for 6*13 obtained from 10^ ran- 
domly generated real and imaginary three dimensional sym- 
metric matrices (Majorana neutrinos). No constraints given 
by Eqs|3]|6]are applied. Gray vertical band represents exper- 
imental values of ^13, EqjTjat the la level. 
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FIG. 3: Frequency probability for 613 obtained from 10^ ran- 
domly generated real and imaginary three dimensional sym- 
metric matrices (Majorana neutrinos). Constraints by Eqsjsj- 
|6] are applied at the 3a level. Gray vertical band represents 
experimental values of &13, Eq[7]at the 3cr level. 



pure simplicity. A modified versions have been consid- 
ered in many papers, see e.g. P^HTB] . 

Let us note that the see-saw neutrino masses are triv- 
ially not well defined if elements of Mji approach zero in 
relation mj^M^^niD, which can happen as elements of 
Md and Mji are generated in the [—1,1] range. Thus we 
will consider only Dirac and Majorana type of neutrinos. 
To disentangle these two types of neutrinos is of great 
importance for neutrino physics [lOj. If we connect the 
mechanism of generation of their tiny masses with some 
heavy states, this issue starts to be important for high 



energy colliders like e+e [211 



e e 



and LHC 



We can extend this discussion to the case where el- 
ements of the neutrino mass matrix involve complex 
phases (random numbers with modulus in the interval 
[0,1] and a phase in a range [0;27r]). It might seem 
that real and imaginary matrices should produce qual- 
itatively the same results: situation differs only as an 
additional degree of freedom (phase) enters in each entry 
of the mass matrix. However, this is not true, as in fact 
tliey differ. This change does not affect the shapes of the 
eigenvalues distributions in a significant way but other 
observables distributions like elements of the mixing ma- 
trix may change. This is exactly the case for 613. 

Figj2]shows the results of numerical predictions for 6'i3, 
additional constraints given by Eqs|3][6] are not applied. 
We can see that the real and imaginary random matrices 
behave qualitatively differently in vicinity of small 9i3 
and real matrices fit better. That frequency probability 
of random complex matrices tends to zero with 613 — > 
can be understood in the following way. For real matri- 
ces CP complex phase is zero or tt and CP is not violated 
when 9i3 = and therefore there is no direct restriction 
on values of ^13. Reversing the argument, if CP symme- 



try is violated (which is the case for randomly generated 
complex mass matrices) then ^13 must be different from 
zero. We found that frequency probability for the phys- 
ical CP violating phase S which emerges from randomly 
generated mass matrices is constant (so this is a trivial 
plot and we do not show it here. From pure mathematics 
(a probability measure) follows that in such case a prob- 
ability to get S equal to zero or tt is zero, so probability 
distributions should tend to zero for ^13 — 0, what can 
be seen on histograms, see solid line in Fig. [2j If Eqs|3][6] 
are applied, a number of mass matrices which cover the 
experimentally interesting values of ^13 decreases. For 
instance, using la cuts in Eqs 3][6 only 4886 real and 
4686 complex matrices remain out of initial 10*" random 
matrices. Using 3a cuts we are left with 31853 real and 
29141 complex matrices, see Fig|3] Comparing between 
Dirac and Majorana cases more mass matrices in Majo- 
rana case survive the cuts. The main conclusion remains 
however the same as in the case where no cuts Eqs|3][6] 
are applied: there are more matrices which reproduce ex- 
perimental data for real matrices and complex matrices 
do not condense in the experimentally interesting region. 

It is interesting that the last experimental data implies 
the presence of a fourth sterile neutrino pSH^ . Usually 
the so called 3-1-1 and 3-1-2 models of sterile neutrinos 
are considered |35j . We then inspect the mass matri- 
ces of dimensions four and five. Here we focus on an 
element |t/i3p and define an effective angle such that 
sin^ ^13 = |t^i3p, so we can compare the results obtained 
from models with different number of neutrinos. 

In a case of more than three dimensional mass matri- 
ces relations among the elements of the mixing matrix 
and experimentally defined mixing angles are more in- 
volved, therefore we will consider this case without cuts 
driven by Eqs 3][7 and without a discussion of the im- 
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FIG. 4: Shifts in the direction of small values of ^13 with 
increasing number of sterile neutrinos. No constraints given 
by Eqs|3|6|are applied. 



FIG. 5: Probability distributions for \Uii\ elements of the 
neutrino mixing matrix when the constraint Eq[7] for &13 is 
applied. The gray vertical band stands for Scr values which 
follow from global fits [32] . 



pact of additional mixing matrix parameters of four and 
five dimensional mass matrices. We proceed as in the 
previous case and generate random matrices and then 
calculate the mixing matrix. In Fig|4]we show results for 
Majorana (symmetric) real matrices of different dimen- 
sions. We observe that the more neutrinos we have the 
bigger number of random matrices reproduce the exper- 
imental value of the ^13 parameter. Similar results give 
Dirac neutrinos. 

We have no clear explanation why the probability of 
getting small values for ^13 grows with dimensionality of 
matrices. There is a mathematical law of Wigner [55] 
connected with real symmetric random matrices which 
states that their eigenvalues accumulate around zero 
(Wigner's semicircle law [30] is exact in the limit of infi- 
nite dimensions of matrices). To our knowledge, there is 
no relation connecting the distributions of the values of 
elements of eigenvectors with the dimensionality of real 
symmetric matrices. 

Nonetheless, let us try to understand the result we ob- 
tained. In the anarchical model of neutrino masses, when 
no cuts on parameters are applied, mixing angles have 
very similar probability distributions. This is because 
every element of the mass matrix has exactly the same 
probability distribution. Therefore, if the dimension of 
the mass matrix grows, the unitarity of the mixing ma- 
trix constrains sums of squared modulus of mixing matrix 
elements, and an average value of each of elements of the 
mixing matrix will decrease with increasing dimensional- 
ity. This forces the probability distribution to take higher 
values for small angles as the dimension grows, a kind of 
pattern obtained in Fig(4] 

We can make a step further and consider for the first 
time implications of nonzero ^13 on values of elements 
|J7i4| and \U24\ in the 3-1-1 sterile model when the con- 
straint Eq[7]for ^13 is applied. (Without cuts, probabil- 
ity distributions for |?7i4| and IC/24I are practically the 
same as for in Figj2| as the just discussed equal- 
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FIG. 6: Probability distributions for \U24\ elements of the 
neutrino mixing matrix when the constraint Eq[7] for &13 is 
applied. The gray vertical band stands for upper bounds de- 
rived by MINOS [33]. 



distribution argument suggests.) It is interesting that 
analyses of recent experiments predict non-zero values 
for them [3T1[35]. In Fig. |5] and Fig. [6] the plots are given 
for real and complex neutrino Majorana mass matrices 
(Dirac patterns are very similar). We can see that proba- 
bility distributions for both real and complex cases reach 
a maximum at larger values of |J7i4| (|J7i4|maa; — 0.7 
for a real case, |?7i4|max — 0.58 for a complex case) 
when compared with \U24\ {\U2i\max — 0.1, real case; 
|t^i4|maa; — 0.4, complcx casc). This is in agreement 
with fits discussed in [23 (|t^/x4p < I^e4p for both low 
and high energy neutrino experimental fits) . In Figjs] we 
can see that the most frequent values of |J7i4| do not co- 
incide with its experimentally preferable region, though 
in this experimental region they can also be substantial. 

More involved analyses including other experimental 
constraints and various mass hierarchies in 3 -I- 1 and 
3 + 2 models go beyond this Brief Report and are left for 
a future work. 

In conclusion, though random matrices can not solve 
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fundamental problems in neutrino physics, they generate 
intriguing hints on the nature of neutrino mass matrices. 
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